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Abstract of EP0749014 

The method involves reconstructing a field, such 
as that of the velocities of particles in a fluid, 
using prior knowledge of the physical phenomena 
in order to produce a regularising model and an 
analytic function. A field network which adapts 
itself to the data available is then created using a 
Kohonen type algorithm. This produces a good 
quality of uniformity tfiroughoutthe network, ^in^ • . 
particular the angles of adjacent sides of a net* 
are greater than a threshold value and there is a 
minimum distance between nodes. The field is 
then reconstructed by optimising the analytic 
function. 
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(57) La pr§sente invention conceme un proced6 de 
reconstruction d'un cliamp de caracteristiques asso- 
ciees k un ph6nomdne physique k partir de mesures 
^parses comportant les Stapes sulvantes : 

elaboration d'un moddle de reconstruction du 
champ elabore a partir de connaissances a priori 
sur le ph^nom^ne dtudi^ de fagon a obtenir une 



DONNEES EPARSES FOURNIES PAR 
LE PROCEDE DE TRAJECTOGRAPHIE 
D'OBJETS 




fonction analytique ; 

realisation d'un maillage du champ ; 

reconstruction du champ par optimisation de la 
fonction analytique ; 

et dans lequel le module de reconstruction prend en 
compte les caracteristiques physiques du phenomene. 
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Description 

Domaine technique 

Le domaine technique auquel se rattache I'lnvention est celui de la reconstruction d'un champ dense de primitives, 
ou caract6ristiques de trajectoires d'objet, comme la Vitesse ou l'acc6l6ration par exemple, k partir de mesures ^parses. 
Plus particulierement, I'invention concerne un proc6d6 appllqu^ au domaine de la v6loclm6trie par imagerie de partl- 
cules en mdcanlques des fluides. 

Etat de la technique ant6rieure 

Induire et reconstruire un champ de primitives 6quiparties k partir de la connaissance de donnees tres ^parses 
est un probleme commun k de nombreuses applications de domalnes varies. II concerne autant la construction de 
cartes topographiques k partir de donn6es issues de relev6s terrains, que la reconstruction d'un champ de vitesses 
dense k partir des resultats fournis par I'^quation de contrainte du mouvement en estimation de mouvement ou encore 
la reconstruction de la surface visible d'objets a partir de mesures st6r6oscopiques. 

En m^canique des fluides. il est important de disposer de la connaissance du champ de Vitesse en tout point du 
domaine 6tudi§. 

Trajectographie d'objets 

Un article de C. Peterson (reference [1]) pr6sente un algorithme neuromim6tique pour trouver tes trajectoires de 
particules dans des exp6riences de physique des particules. L'algorlthme est bas6 sur un r6seau de neurones du type 
de celui de Hopfield et sur des equations de mise k jour des neurones issues de la th^orie du champ moyen. Le 
probleme de suivi de particules est formalise comme la recherche de trajectoires (objets) compos6es du plus grand 
nombre de points (fragments) formant des courbes lisses et droites. Un neurone du rSseau 6tant associe a chaque 
couple de fragments, ie probleme revient k rechercher des courbes formdes de couples de fragments contigus, de 
m§me sens et le plus alignes possible. Une 6nergie E, fonction pond6r6e des temnes assocles k ces crlt^res codant 
la qualite de la solution, permet d'obtenir les equations d'evolution des coefficients synaptiques du reseau neuronal. 
A la convergence du reseau, Tenergie E est optimisee. On peut noter les defauts suivants : 

une grande sensibilite aux parametres. Le resultat final est tres dependant des paramdtres pond6rant les diffdrents 
termes de I'^quatlon de I'^nergie ; 

non satisfaction des contraintes. Le reseau optimisant une fonction, qui code globalement le probleme, on ne peut 
pas assurer que toutes les contraintes seront rigoureusement satisfaites. En d'autre termes, la solution peut con- 
tenir des elements ne satisfalsant pas les contraintes. La solution finale peut done etre de mauvalse quality ; 
temps de calcul prohibitif. L'algorlthme de calcul ne realise aucun pretraitement. Si I'image contient n attribute, le 
reseau sera compose de rfi neurones, soit poids synaptiques k determiner. Les temps de calcul requis sont 
done tres importants. C'est d'ailleurs la raison pour laquelle cet algorithme n*a 616 implants que sur des problemes 
de tres petite dimension. 

Un article de A. Cenedese. G. Romano, A. Paglialunga et M. Terlizzi (r6f6rence [2]) presente une methode neu- 
romim6tique de reconnaissance de trajectoires en m6canique des fluides. Le r6seau utilise est un reseau multicouches 
base sur un module fourni par le classifieur de Carpenter-Grossberg. Dans ce type de reseau. le calcul des poids 
synaptiques s'effectue sans supervlseur. L'algorlthme fonctionne avec une sequence d'images. L'apprentissage. et 
done le calcul des poids synaptiques, est realise grace k un certain nombre d'images de la sequence. Ensuite le champ 
de Vitesse pourra etre determine, Le reseau neuromlmetique d6crit dans cet article est un reseau k apprentissage. 
L'inconv6nient majeur de cette mdthode est done que, pendant le temps d'enregistrement de la sequence, le regime 
du fluide en mouvement analyst doit etre stationnaire. D'autre part, la methode necessite un temps de calcul trks 
Important et utilise de nombreux parametres dont I'estimation est ad hoc. Enfin, en l'6tat actuel, l'algorlthme propose 
necessite une importante validation. 

Un article de M. Gyulassy et de M. Harlander (r6f6rence [3]) propose un algorithme de suivi elastique pour deter- 
miner des trajectoires d'exp6riences de physique des particules. Dans I'approche de suivi 6lastique, une trajeetoire 
est un objet h6licoidal, qui se trouve dans une forme s'ajustant fe mieux aux donnees, i.e. aux fragments. L'helice peut 
§tre vue comme etant chargee 6leetriquement et 6XanX attlree par les fragments qui ont une charge de signe oppose. 
La notion de suivi 6lastique intervient dans le calcul des poids synaptiques. Le r6seau de neurones utilise est le meme 
reseau que celui utilise par C. Peterson dans I'artiele precedemment cite. Seul le calcul des poids synaptiques est 
modifie, chaque poids synaptique tenant alors compte de la qualite de la trajeetoire ajustable sur les neurones relies 
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par ce poids synaptique. Cette mdthode pr6sente les mdmes lnconv6nients que la mdthode de C. Peterson : 

sensibility aux paranri^tres, 
non satisfaction des contraintes, 
5 - temps de calcul prohibitif. 

La trajectographie d'objets conceme I'^tude des trajectoires tennporelles et spatiales d'objets en nnouvement, par 
exemple des particules en nnouvement. Elle identifie et reconstitue les trajectoires des objets et fournit les caractdris- 
tiques de ces trajectoires (par exemple vitesse, acceleration). Comme les objets se d6placent dans un certain domains 
10 d'etude et comme un proc^de de trajectographie d'objets ne fournit les caract^ristiques des trajectoires que de ces 
objets, les caract^ristiques des trajectoires ne sont connues qu'en quelques positions du domaine etudie. Ces carac- 
teristiques de trajectoires. que Ton appelle mesures ou donn^es, sont done ^parses. Le r^sultat du proc^de de trajec- 
tographie d'objets va constituer I'entrde du proc6d6 de reconstruction d'un champ dense de primitives, i.e. caractdris- 
tiques de trajectoires d'objets. 

^5 La reconstruction d'un champ dense de caract6ristiques de trajectoires d'objets concerne la reconstruction d'un 

champ dense oij les caracteristiques des trajectoires d'objets sont connues en tout point du domaine 6tudl6 h partir 
des mesures ^parses fournies par le precede de trajectographie d'objets. 

Les caracteristiques des trajectoires d'objets fournies par le procede de trajectographie d'objets doivent 3tre com- 
bin6es et dlffus6es pour obtenir la valeur des caracteristiques de trajectoires en tout point du domaine etudie. Ce 

20 probieme s'appelle la reconstruction d'un champ dense de primitives. II ne s'agit pas d'un problems d'interpolation de 
donnees, mais d'un problems d'approximation. En effet, les mesures ne peuvent dtre consider6es comme absolument 
fiables. Des erreurs au cours du traitement preaiable ont pu apparattre. Du bruit peut par ailleurs alterer les resultats. 
II ne faut done pas vouloir interpoler les mesures connues exactement: contraindre les donnees k conserver leurs 
vaieurs conduirait inevitablement a un resultat errone. La condition d'interpolation doit dtre relaxee, ce qui transforms 

25 le problems en un probieme d'approximation. . 

Methodes de reconstruction utillsees en mecanlque des fluldes 

Des methodes ayant pour objet la reconstruction d'un champ dense de caracteristiques de trajectoires sontdecrites 
30 dans differents articles de I'art connu en mecanique des fluides, et notamment dans ceux cites ci-apres. 

Une methode couramment utilisee (reference [4]) utilise la methode d'interpolation de donnees developpee par 
P. Hardy (reference [5]) pour des applications en topographie geophysique. L'hypothese falte est que le champ de 
Vitesse k reconstruire peut etre assimiie k une surface multiquadrique : 



35 
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K =(''M^i) au point (x.,^,) 



45 ou Vj est le vecteur vitesse au point (X/, yj) , Uj et v,- sont les composantes du vecteur V/,A/ est le nombre de vecteurs 
Vitesse connus et aj et bj sont des constantes a determiner En appliquant ces relations sur chaque vecteur vitesse 
connu, on obtient equations k 2A/inconnues (coefficients aySt bp. Les constantes inoonnues ayet bjsonx detemninees 
en resolvant le systeme d'equations lineaires trouve. Les constantes ainsi calcuiees sont alors utillsees conjointement 
avec les vaieurs Xy et associees, pour reconstruire le champ de vitesse en tout point {x,,yD de I'ecoulement, gr§ce 

50 aux equations (1 ). L'Inconvenient majeur de cette methode reside dans le choix arbltraire de relations multiquadriques : 
aucune etude n'a ete realises permettant de justifier theoriquement ou experimentalement le choix de telles relations 
dans le cadre d'ecoulements de mecanique des fluides. 

Agui et Jimenez (reference [6]), apres avoir compare diverses methodes d'interpolation incluant Tinterpolatlon 
polynomiale, les moindres carres polynomiaux et les convolutions, ont propose une methode de convolution simple 

55 avec une fenetre Gausienne adaptative. La vitesse au point (x^ yj) est ainsi deflnie par: 
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u,{resp.v^) = 



Zj=, aj.Uj(resp.v,) 



Les coefficients ay sont des coefficients d' Interpolation d6finis par: 



a J = exp 



La largeur de fenetre optimale H propos6e par Agui et Jimenez est H= 1 ,24.5. ou 5 est la distance moyenne antra 
particules calcul6es par 5 = (l/^cn)^ avec n nombra moyen de particules par unit6 de surface. Comma dans la cas de 
la precedente methode, rien ne prouve qu'interpoler par convolutions a una quelconque signification physique par 
rapport aux propri6t6s de i'^coulement de mecanique des fluides. En consequence, s'il est vrai que le champ de Vitesse 
reconstruit pr6sentera un aspect qualitatlvement "esthetique", grace au llssage sous-jacent, 11 n'en demeure pas molns 
qu'll n'aura aucun sens physique r6el et qu'll sera quantltativement Inexploltable. 

Una interpolation lineaira par methode das moindras carr^s a et6 proposea par K. Iwaichi at K. Ohmi (reference 
[7]). En supposant qu'll existe deux equations lineaires entre tes vitesses inconnues Vf^u^vj^ et les vitesses connues 
VJ(t/y,i^y)ouyG<1.A/>: 



(2) 



Les variables Inconnues (u-pViBX leurs d6riv6es) sont evalu6es en tout point (X/,y3 k partlr des obsen^ations connues 
^^j'^j\^jyj> ^® sorte que la difference reliant les deux termes des egalltes (2) soft mlnimis6e. Une mfime cause 
ayant une meme cons6quence, cette Interpolation 6tant Iin6aire. elle fait abstraction de toutes les proprl6t6s physiques 
des ^couiements. 

Enfln. HIno et al. (reference [8]) ont expose une methode pour estimer le champ de vltesse Instantan6 k partlr de 
quelques points mesur6s. Dans une premiere 6tape, la methode permet d'interpoler les donnees connues en consl- 
d^rant que le champ de vitesses k reconstruire est la deformation d'une hypoth§tique plaque mince 6lastique par 
quelques charges virtualles concentrees en des points arbitraires (la deformation de la plaque correspond aux com- 
posantes de la Vitesse uet v). L'utilisation de plaques minces ou membranes est largement repandue en regularisation. 
La methode Hino et al. a cependant la particularity d'utiliser le fomr^alisme de la fonctlon de Green associ^e k la 
deformation de la plaque mince pour resoudre le probleme. Dans une seconde etape, afin d'am6Morer Testimation 
grossiere foumie par la premiere etape, la fonctlonnelle sulvante est minimisee : 



ou Q est le domalne de l'6coulement, Uq etVo sont les estimations de u et \/fournies par la premiere etape et X est le 
multlplicateur de Lagrange. Le premier terme a pour objet de s'approcher de restlmation grossiere dej^ calcuiee, tandis 
que le second introduft Cequation de continuite de m6canique des fluides (div(V)=0). Par rapport aux autres m6thodes, 
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cette m^thode realise une approximation et non une interpolation et s'affranchit done ainsi du probl§me \\6 au bruit ou 
aux erreurs. D'autre part, elle introduit une propri6t6 de m6canique des fiuides.. 

Matgrd des avancdes marquantes, les m^thodes d'approximation de donndes sont encore loin de foumir des 
r6sultats satisfaisants, notamment en ce qui concerne les applications en reconstruction de champ de Vitesse dense. 
s Si Ton cherche ^ faire le point sur les carences des m^thodes traditionnelles, on peut se rendre compte que : 

les mdthodes utilisdes, principalement en m^canique des flu Ides, ne prennent que rarement voire jamais en compte 
des considerations sur le ph^nom^ne physique sous-jacent ; 
la plupart du temps, I'a priori introduit est bas6 sur des termes de lissage ; 
10 . revaluation de la quality d'une reconstruction est blen difficile k r§aliser; 
la determination des param^tres reste une question ouverte ; 
les temps de calcul requis par ces m6thodes sont en general tres importants ; 
aucune de ces mdthodes ne prend en compte la repartition spatiale des mesures. 

Ces procedes de I'art anterieur ne permettent pas d'obtenir de reconstruction de champ dense de caracteristiques 
de trajectoires d'objets ayant un reel sens physique. 

LMnvention a pour objet une nouvelle formalrsation qui prend en compte les proprietes physiques d'un phenomene 
aboutissant k un probleme dit d'optimisation, et resolution du probieme sur un maillage de champ optimise et tenant 
compte de la repartition spatiale des objets. 
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Expose de rinvention 



La presente invention propose un procede de reconstruction d'un champ dense de caracteristiques (par exemple 
de trajectoires d'objets, comme la vitesse ou I'acceieration) associees k un phenomene physique a partir de mesures 
25 eparses, caract6rise en ce qu'il comporte les etapes suivantes : 

elaboration d'un modeie de reconstruction du champ eiabore k partir des connaissances a priori sur le phenomene 
etudie de fa9on k obtenir une fonction analytique (appeiee energie k optlmiser) ; 
realisation d'un maillage du champ ; 
30 - reconstruction du champ par optimisation de la fonction analytique ; 

et en ce que le modeie de reconstruction prend en compte les caracteristiques physiques du phenomene. 

L'etape d'eiaboratlon d'un modeie de regularisation formalise les specificites du probieme de reconstruction. Ce 
modeie est elabor6 a parti de connaissances a priori sur le phenomene etudie. It conduit k la definition d'une fonction 
3S analytique, appeiee energie, k optimiser. 

Dans retape de maillage, le maillage est un maillage auto-adaptatif du domaine etudie, adapte k la repartition 
spatiale des mesures foumies par un procede de trajectographie et reparties non-uniformement dans le domaine 
etudie. Cette etape foumit une grille composee de noeuds et de mailles. 

Dans cette seconde 6tape on utilise un algorithme de type Kohonen dans lequel, afin d'assurer une bonne qualite 
40 de regularisation sur le maillage, les mailles de celui-ci dont aussi regulieres que possible, et en ce que des contraintes 
sont imposees sur les angles entre ressorts d'une nnaille : les quatre angles formes par les ressorts pris deux par deux 
devant etre superieurs k un angle minimal. 

De plus, on impose une distance minimale inter-noeuds. 

Dans une realisation avantageuse les caracteristiques sont les caracteristiques de trajectoire d'objet telles que la 
45 Vitesse ou I'acceieration. Les mesures eparses sont fournies par un precede de trajectographie d'objets. 
Le procede de trajectographie d'objet comprend les etapes suivantes : 

une etape d'enregistrement de signaux composes de fragments caracteristiques des positions des objets k diffe- 
rents instants, de parties qui presentent la meme nature que les fragments mais qui sont du bruit et de parties 
50 attribuables au bruit sans confusion possible ; 

une etape d'extraction de parties de signaux ayant la mSme nature que les fragments et determination de carac- 
teristiques et leur etant associees, k partir des connaissances a priori sur leur nature et sur les trajectoires 
recherchees ; 

une etape de repartition de I'ensemble des parties de signaux precedemment extraits en classes, chaque classe 
55 representant une trajectoire potentielle et comprenant un nombre predetermine de fragments ; 

une etape de selection d'un sous-ensemble de classes qui satisfait des contraintes liees au type de reseau obsen/e, 
comportant : 
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• la g6n6ration d'un ensemble de n-uplets constituds d'un sous-ensemble de n classes, 

• la mesure de compatibility des classes assemblies en n-uplets par une premiere fonction analytique, d6ter- 
mrn6e k partir des contraintes lides au type de ph6nomene observd et au dispositif d'acquisition, 

• la mesure de la quality des n-uplets par une second fonction analytique, determinde k partir des contraintes 
liees au phinomene observe et au dispositif d'acquisition ; 

une etape de selection parmi les classes representant les trajectoires potentielles de celles satisfaisant les con- 
traintes, par un procddd d'optimisation sous contrainte du (ou des) crlt6re(s) precedent(s). en utilisant les mesures 
pr6c6dentes de compatibility et de quality, de fa9on k obtenir les trajectoires "ryelles" des objets. 

Le procedy de I'invention est applicable dans le domaine de la vyiocimytrie par imagerie de particules en myca- 
nique des fluides. 

La fonction analytique peut alors etre telle que : 

£(/)=Z,^.-)||/'*'-/t 

ou : 

• ^ 6t designent respectivement les mesures de Vitesse aux sites s et / du champ, 

• et ^^dysignent respectivement les primitives (valeurs de vitesse) k reconstruire aux sites s et f du champ, 

• \i est la viscosity dynamlque du fluide considere, 

• Xj., X\\, XV et p{dg^ sont des paramytres, 

• D est Tensembfe des mesures, 

• S est I'ensemble des sites du champ, 

• n^f est le vecteur unltaire dirigy de s vers f, est un vecteur unitaire normal k et d^^ est la distance syparant 
les sites s et U 

• V.r est la divergence de la fonction f . 

Une fois la fonction analytique E(/>definie, la solution f6u probleme est obtenue lorsque E(f) est optimisye, i.e.: 

^ = argument optimum de f^^J^if) 

ou Fest I'ensemble des solutions possibles. 

L'invention propose un procydy plus robuste car il prend en compte ; 

la nature du phynomine ytudiy ; 

la rypartltion spatiale des mesures utillsyes comme entrye du procydy. 

Par la formallsation du problime sous forme d'un probleme de regularisation, par ryiaboration d'un mod6le dyfini 
k partir de connaissances a priori sur le phenomene etudiy et par I'emploi d'un maillage adapty k ta rypartition spatiale 
des mesures, on peut garantir une bonne quality de la solution obtenue, quality jusqu'^ prysent non atteinte par les 
autres mythodes de t'art connu. 

Par rapport k Tart connu, le precede de l'invention reunit les avantages suivants : 

il utilise les propriytys et caractyristiques du phynomene ytudiy. La solution est done dirigye de fagon juste. La 
solution possyde une ryelle signification au sens de la physique du phynomfene ; 

le maillage auto-adaptatif s'ajustant k la repartition spatiale des donnies permet de stabiliser numeriquement le 
probleme (en termes mathymatiques, on parle d'ameliorer le condition nement numyrique du systdme lie au 
probiyme) ; 

le maillage auto-adaptatif permet d'accyiyrer considyrablement les temps de calcul. 
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Le proc^d^ de i'invention est particuli^rement applicable dans le domaine de v6locim6trie, par exemple en m6ca- 
nique des fluides. Dans cette application particuii^re les modules de I'art antSrieur ne prennent pas en compte la 
physique du ph^nomdne ^ cause de la complexity des calcuts. L'homme de I'art de la vdlocimdtrie n'est pas assez 
math^maticien pour connaltre suffisamment d'outils permettant de mieux colier k la r^alit^ des choses. Jusqu'^ present 

5 ce type de r^sultats lui suffisait. Cela n'est plus le cas actuellement, d'ou t'int^rdt de I'invention. 

Le maillage adaptatif est utitisS dans les probldmes d'analyse de donn^es. L'utllisation d'un maillage adapts n'est 
pas applique dans le domaine de la v^locimdtrie parce que les gens n'avaient pas I'utilitd de rdsultats de meilleure 
quality, nl dans celui de I'optimisation de fonctions parce que les probldmes k rdsoudre ne sV pretaient pas (pas de 
repartition spatiale des mesures). 

10 Le probl6me r6solu par I'invention est celui d'adapter un algorithme d'optimisation au maillage adaptatif. Les dlf- 

ficult6s rencontr6es sont par exemple : des modifications de param^tres tels que la prise en compte de la surface des 
mailles, la distance inter-noeuds dans le cas de la veiocimetrie et dans le cas g^n^ral la mise k jour des sites du champ. 
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Br6ve description des desslns 

La figure 1 illustre le synoptique du precede de invention ; 

la figure 2 illustre ('architecture d'une carte topologique du Kohonen ; 

la figure 3 illustre des representations des forces tangentielles et normales ; 

la figure 4 illustre des representations des projections du vecteur Vitesse Vs \ e\\ et e±. 

Expose detailie de I'invention 



Dans la suite de la description nous considererons k titre d'exemple le domaine de la trajectographle d'objets. 
Le precede de I'invention constitue alors la suite d'un procede d'obtention de trajectoires d'objets en mouvement, 
25 optimisant au moins un crit6re issu de la physique du phenomdne observe, qui peut comprendre les etapes suivantes : 

une etape d'enregistrement de signaux composes de fragments caracteristiques de positions des objets k diffe- 
rents instants, de parties qui presentent la mdme nature que les fragments mais qui sont du bruit et de parties 
attribuables au bruit sans confusion possible ; 
30 - une etape d'extraction de parties de signaux ayant la meme nature que les fragments et determination de carac- 
teristiques leur etant associees, k partir des connaissances a priori sur leur nature et sur les trajectoires 
recherchees ; 

une etape de repartition de I'ensemble des parties de signaux precedemment extraits en classes, chaque classe 
repr6sentant une trajectoire potentielle et comprenant un nombre predetermine de fragments ; 
^ - une etape de selection d'un sous-ensemble de classes qui satisfait des contraintes liees au type de reseau observe, 
comportant : 

• la generation d*un ensemble de n-uplets constitues d'un sous-ensemble de n classes. 

• la mesure de la compatibilite des classes assembiees an n<uplets par une premiere fonction analytique, de- 
^0 terminee k partir des contraintes liees au type de phenomene obsen/e et au dispositif d'acquisition, 

• la mesure de la qualite des n-uplets par une seconde fonction anatytique, determlnee k partir des contraintes 
liees au phenomene observe et au dispositif d'acquisition ; 

une etape de selection parmi les classes representant les trajectoires potentielles de celles satisfaisant les con- 
45 traintes, par un procede d'optimisation sous contraintes du (ou des) critere(s) precedent(s). en utilisant les mesures 

precedentes de compatibilite et de qualite, de fagon k obtenir les trajectoires "reelles" des objets. 

L'entree du procede de {'invention est constituee par la sortie de ce procede d'obtention de trajectoires d'objets 
en mouvement. Le procede de I'invention comprend trols etapes, illustrees sur la figure 1 : 

60 

dans une premiere etape 10, un modeie de r6gularisation est eiabor6. Celui-ci prend en compte des caracteristi- 
ques et proprietes sur le phenomene etudie. Une fonction analytique E(f) est definie. La solution f du probieme 
est obtenue lorsque cette fonction analytique E(f) est optimisee ; 

une seconde etape 11 consiste k realiser un maillage du domaine etudie. maillage qui est adapte k la repartition 
^ spatiale des mesures ; 

une troisieme etape 12 procede k Toptimisatlon de la fonction analytique E(f) definie en etape 10 sur la grille 
definie par le maillage realise en etape 11 . 
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L'invention consists k combiner les differentes m^thodes ou parties de m6thodes existantes dans I'arl connu et 
les adapter pour les appliquer k un probleme particulier, h savoir la reconstruction d'un champ dense de primitives k 
partir de mesures ^parses. 

L'interet du procede de l'invention par rapport k Van anterieur est de proposer une m6thode globale qui, pour plus 
complexe qu'elle soit, fournit des resultats plus precis, plus robustes au bruit sur les mesures et possedant une r6elle 
signification physique au sens du phenomene. 

On va h present etudier chacune des stapes du proc6d6 de l'invention. 

Etape d'^laboratlon du mcxidle de regular Isation 

Le probleme d'approximation de donn6es, comme la grande majorlte des probl6mes rencontres en traitement 
d'images, correspond k un probleme inverse : a partir des observations du comportement d'un systeme et a partir 
d'informations a priori, on cherche k remonter k une information non directement observable et caracterisant ce sys- 
teme (le champ dense de vitesses instantan6es). On parle de probleme inverse, puisqu'il s'agit d'inverser les Equations 
decrivant le comportement du systeme. Or, si I'etablissement des equations directes ne souleve pas d'autres difficultes 
que celles de ia mod6lisation de processus physiques plus ou moins complexes, il n'en est pas de m§me pour {'inversion 
de ces m§mes Equations. C'est ainsi que beaucoup de ces probl6mes inverses appartiennent k la cat6gorie des pro- 
blames "mal pos6s". C'est k Hadamard que Ton dolt la notion de "problfemes correctement pos6s", introdulte dans le 
but de savoir quels types de conditions aux limites sont les plus naturels pour diff6rents types d'6quations diff 6rentielles. 
Ainsi est dit probleme mal pos6, un probleme qui ne verifie pas une des trois conditions donn^es par Hadamard pour 
qu'un probleme soit dit bien pos6, a savoir qu'une solution existe pour toute observation, que cette solution soit unique 
et enfin qu'elle depende continument des donn^es. 

Le probl6me de reconstruction d'un champ dense de primitives est un problfeme mal pose au sens de Hadamard. 
En effet, ce probleme est caracterise par l'incompl6tude des observations. Les donn^es incompletes r^parties non- 
uniformement sur la surface k analyser restreignent le champ de vitesse localement, mais ne le determinent pas de 
fagon unique en tout point de la surface, it existe done une infinite de champs de vitesses possibles. La th6orie de la 
regularisation vise precisement ^ transformer le probleme mal pose en un probleme bien pos§ en restreignant la classe 
des solutions admissibles par I'lntroduction de connaissances a priori sur celles-ci. En effet, en I'absence de conside- 
rations compl6mentaires, la seule mod^lisatlon du lien entre les observations du comportement du systeme et les 
caracteristiques de la solution a trouver ne permet pas de reconstruire une approximation tant soit peu satisfaisante 
de la solution. Pour essayer de contoumer cette difflculte, il est n^cessaire d'emettre des hypotheses sur les propri6t6s 
de la solution et de les int6grer dans le processus de reconstruction. Cette introduction de contraintes suppl6mentaires 
liees a une connaissance a priori est qualifiee de regularisation. 

La resolution du probldme de reconstruction d'un champ dense de primitives passe done par la regularisation, i. 
e. par la recherche d'une solution approchee, unique et stabilisee par I'lntroduction explicite d'une information a priori. 
Des renseignements sur la theorie de la regularisation peuvent etre trouves dans les references [8] et [9]. 

L'information a priori peut etre de nature deterministe ou stochastique, ce qui permet de distinguer deux types 
d'approches : 

techniques de regularisation deterministes : voir references [11] et [12] ; 
techniques de regularisation stochastiques: voir references [13] et [14]. 

L'une ou I'autre de ces techniques revient k d6finir une fonction analytique k optimiser et k formaltser le probl6me 
comme un probleme d'optimisation k resoudre. La fonction analytique a optimiser est toujours constituee de deux 
termes. Le premier, appele energie d'attache aux donnees, mesure I'adequation entre les mesures et le champ de 
primitives reconstruit. Le second, appeie energie contextuelle, decrit toutes les connaissances a priori sur le champ 
de primitives k reconstruire. 

La fonction analytique E(f)k Mfinir peut contenir des connaissances liees au processus de formation des images, 
au bruit envisage sur les mesures et surtout sur la nature du phenomene etudie. Les connaissances a priori sur la 
nature du phenomene etudie pourront etre par exemple : 

des contraintes de regularite: en tout point du domaine etudie, les variations de valeurs de primitives doivent etre 
faibles ; 

des contraintes liees k I'equation de contrainte de mouvement en mecanlque des fluides. Si la primitive k recons- 
truire est le vecteur vitesse, on peut imposerque divVsoiX nul en tout point du domaine etudie, ou d/u'V^represente 
la divergence de la vitesse V ; 

des contraintes liees aux forces de pression et de frottement dans le cas ou la primitive a reconstruire est le vecteur 
Vitesse dans un ecoulement de mecanique de fluides. 
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• 



Dans le cadre de rapplication particuli^re de i' invention au domaine de la v6locim§trie par imagerie de particules 
en mdcanique des fluides, la fonction d'energie cholsle est : 

ou : 

IS • et d6signent respectivement les mesures de Vitesse aux sites s et fdu champ, 

• et ff d^signent respectivennent les primitives (valeurs de vitesse) a reconstruire aux sites s et f du champ. 

• |i est la viscosity dynamique du fluide consld^rd, 

• XJ_, A.\\, V7 et p{dgf) sont des param6tres, 

• D est Tensemble des mesures, 

20 • s est I'ensemble des sites du champ. 

• n^f est le vecteur unltaire dirig^ de s vers t, est un vecteur unitaire normal k n^f et d^f est la distance s^parant 
les sites s et t, 

• V. / est la divergence de la fonction f . 

25 En annexe, on trouve une explication justifiant le choix de cette expression pour E(/). 

Une fois cette fonction analytique E(/) d6f inie, la solution f du probl^me est obtenue lorsque E(f) est optlmiste, i.e. : 



? = argument optimum de y^pE(f) 

30 

ou F est I'ensemble des solutions possibles. 
Etape de malllage adaptatif 

3S Les approches de regularisation conventionnelles commencent par I'^laboration d'un module d'energie en r6alisant 

souvent une formulation variationnelle du probldme. Puis celle~ci est suivie par une discretisation de la fonctionnelle 
d'dnergie et des contraintes aux limites sur une grille de calcul r^gulidre couvrant le domaine dtudid. 

Dans le cadre de la reconstruction k partir de mesures ^parses k laquelle a trait le present procddd, il apparaTt 
opportun de prendre en compte la repartition statistique spatiale des donnSes sur le domaine 6tudi6 k {'aide d'un 

^0 maillage adaptatif. 

Dans les zones du domaine dtudi^ ou I'on n'a pas ou peu de donndes, une grande precision n'est pas n^cessaire, 
dans la mesure ou seul le terme d'6nergie contextuelle sera pris en compte. A I'inverse, il est important d'avoir beaucoup 
de mailles aux endroits ou la density de donndes est forte. A la place d'un malllage rdgulier, il est done ndcessaire 
d'utiliser un maillage adapts k la density spatiale des donnees. 
45 Un autre avantage reside dans le gain en temps de calcul. En effet, utiliser des maillages k mailles uniform6ment 

r6parties peut conduire k une tr6s grande lourdeur calculatoire, parce que, afin d'avoir un nombre de mailles assez 
grand dans les zones aux variations importantes, il est n^cessaire d'utiliser un maillage tr6s fin sur I'ensemble de 
rimage et notamment aux endroits oD cette finesse ne serait pas ndcessaire. 

Enfin, on a montr6 que I'emploi d'une grille non r6guli6re qui s'ajustait aux donn6es. de telle sorts que chaque 
so maille de la grille tende k contenir au moins une donnee, permettait d'am6liorer le conditlonnement num6rique du 
syst^me Ii6 k Toptimisation de E(f) et done d'am6liorer la stabllit6 de la solution par rapport aux bruits sur les mesures. 

On va k present considerer la m^thode de maillage adaptatif. Kohonen (voir reference [1 5]) a d^veloppe un module 
de r^seau de neurones s'appuyant sur {'observation neurophysiologique d'autOK>rganisation des r^eaux de neurones 
biologiques rdels. Ce r^seau de neurones est capable de s'auto-organiser de telle sorte qu'il approxime la fonction de 
ss density des vecteurs d'entr6e du r6seau de neurones. 

Ce type d'algorithme d'auto-organisation est utilise pour des applications varices : reconnaissance de la parole 
(voir r6f6rence [16]), traitement du signal (voir r6f6rence [1 7]), probi^me de th6orie des grappes (voir rdfdrences [18], 
[19] et [20]), applications industrielles (voir rdfdrence [21], choix de la dimension d'une representation rdduite d'un 
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ensemble de points de R". construction de representations hi6rarchis6es. Des applications electronlques de la carte 
de Kohonen ont 6te proposees (voir references [22] et [16]). 

Par rapport k I'algorithme classique, plusleurs modifications ont dO §tre apportees au reseau de neurones de 
Kohonen. 

Dans I'algorithme classique de Kohonen, on suppose que I'espace des stimuli (entrees) est structure, qu'on peut 
definir une notion de voisinage ou de ressemblance, par exemple au sens des voisinages d6finis k partir de la distance 
euclidienne de R". On suppose que I'on dispose de N unites (ou neurones de sortie), sur lesquelles est defini un 
syst^me de voisinage et on note V(i) le voisinage de Vun\X6 i (y compris I'unite /). On suppose plus precisement que 
I'espace T des stimuli est contenu dans I'espace euclidien R". 

Le r6seau de neurones defini par Kohonen est constitue d'une couche de neurones d'entr6e (analogiques ou 
binaires) d6finissant I'espace des stimuli possibles et d'une couche de neurones de sortie mise en correspondance 
avec I'espace des prototypes (voir figure 2). 

La rdponse de I'unlt6 /de sortie k la presentation d'un stimulus v= [v^, .vj^est donn6e par : 



ou les eR sont des param^tres adaptatifs. En fait, chaque ensemble ordonne (^^,..,.,jij peut etre regard^ comme 

une sorte d'image pouvant etre comparee ou associee k I'ensemble ordonne correspondant (v^ v^y Ce qui carac- 

tdrise done I'unit6 /, c'est le vecteur des poids synaptiques #73,-= [^;i,....,fi Jt; repr6sent6 dans le mSme espace Y que 
les stimuli. C'est sur ce vecteur que porte i'apprentissage. 

L'objectif de I'invention est de definir un processus adaptatif d'apprentissage, dans lequel les param^tres de toutes 
les unites convergent vers des valeurs telles que chaque unite devienne sp6cifiquement sensible a un domaine par- 
ticulier des signaux d'entree dans un ordre regulier A cette fin. pour un 6tat donn6 des vecteurs m^, est d6finie Tappll- 
cation 0{m), qui, au stimulus v, associe I'unit6 telle que m^^soit le plus proche de vdans T, ia relation de proximite 
adoptee pouvant dtre la distance euclidienne entre vecteurs. II s'agit 1^ d'un processus de competition: election du 
vecteur de poids le plus proche du stimulus presente. Le probl6me est alors de definir un processus d'adaptation des 
vecteurs m^de sorte qu'ils convergent vers des valeurs m/' qui rendent I'appiication 0(m*) respectueuse des relations 
de voisinage. Pour ce faire, I'apprentissage se deroule de telle sorte que les unites ne soient pas atfectees par I'ap- 
prentissage de maniere independante. Au cours d'une iteration, on d6finit un voisinage v(y autour de I'unitd ga- 
gnante, Seules les unites appartenant a v(y sont affect6es par I'apprentissage. 

T Kohonen a propose I'algorithme suivant (voir r6f6rence [15] : 

1. Au temps t=0, les valeurs de sont quelconques (aleatoires et independantes par exemple). 

2. Au temps t+1 , on presente le stimulus v(^*i), choisi aleatoirement dans T selon une loi de probabilite p{v). On 
determine alors I'unitd de meilleur r^ponse verifiant: 



ce qui est equivalent k : 



||v<-)-ml:)|=min|v(-)-;„(')| 



3. On renforce alors i'ajustement en /^et ses voisins en posant: 

;i^-)=^<'>+a<'>.(v;.-) -;.<')). si i ^V{Q 
= si i^ViK) 
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Le r6sultat 6nonce par Kohonen est que !e vecteur m, 6es coefficients synaptiques converge (au cours du temps) 
vers des vateurs mf*Xe\\es que rapplication <I»(/rf ) respecte la topologie. Et plus pricisdment, la fonction density ponc- 
tuetle des vecteurs m/ converge vers une fonction qui est une image discrdtisde de la densitd p(v) avec laquelle sont 
pr6sent6s les stimuli v, et reflate done ainsi la distribution de probability des entries. 
s Les poids synaptiques reliant les neurones de la couche d'entr^e k un neurone de la couche de sortie d^finissent 

les coordonnees dans I'espace des entries du prototype repr^sent^ par ie neurone de sortie. Apr^s apprentissage, 
un stimulus d'entrde sera reprdsente par ie neurone prototype dont il se rapproche le plus. Par consequent, une fots 
I'apprentissage acheve, les valeurs de connexions ddfinissent un pavage de \A3ronoi de I'espace des entries qui dolt 
^chantiilonner le mieux possible la distribution de probability des stimuli d'entr^e. Ce sera le rdle de I'apprentissage 
10 d'etablir un rapport entre la topologie de i'espace des entries et ia proximity des repr^sentants des prototypes sur la 
couche de sortie. On pense ainsi faire acqudrir au rdseau des capacit^s de robustesse, de stability et de simplicity de 
la reprysentation. 

Dans I'invention on apporte les modifications suivantes : 

*5 • Afin d'assurer une bonne quality de rygularisation sur le maillage, les mailles de celui-ci doivent dtre aussi rygu- 
nyres que possibles (en Toccurrence, les yiements quadrilatyraux doivent ressembler k des carrys autant que 
possible). Les yiyments singuliers (rectangles trop plats, rectangles dygynyrys) pourralent causer des instabilitys. 
Des contraintes sont done imposyes sur les angles entre ressorts d'une maille : les quatre angles formes par les 
ressorts pris deux par deux doivent dtre supyrieurs k un angle minimal. Si, apr^s mise k jour d'un noeud, parmi 

20 les quatre mailles auxquelles il appartient, une maille prysente un aspect plat, la mise k jour de ce neurone n'est 

pas effectuye : le noeud ne change pas de place. 

• II peut ytre nycessaire d'imposer une distance mtnimale inter-noeuds. En consyquence, on impose une seconde 
contrainte sur les distances inter-noeuds. Si. aprys mise a jour d'un noeud, parmi les quatre ressorts auxquels il 
2S est reliy, un est de longueur infyrieure a d^j^, on ne ryalise pas la mise ci jour du neurone concerny. 

Etape de resolution 

L'ytape 1 0 a conduit k la dyfinition d'une fonction analytique E(/) a optimiser, telle qu'a I'optimum de cette fonction 
30 corresponde la solution optimale. 

L'ytape 11 a fourni une grille non ryguli^re adaptye k la density des donnyes. 

L'ytape 12 est une ytape de rysolution au cours de laquelle la fonction analytique E(f) est optimisye sur la grille 
non ryguliy re gyny ry e par l'ytape 1 1 . 

Toute mythode d'optimisation existant dans I'art connu pourra etre utilisye. 
35 Le rysultat fourni par cette ultime ytape est une carte ou les valeurs de la primitive k reconstruire sont connues 

en tout point du domaine ytudie. 

ANNEXE 

40 Justification du choixde {'expression E(f)daM lecas de i*application de veiocimetrieen mecanlque des fluides 

Soit une lattice rectangulaire Sde sites ponctuels correspondant k I'image numyrique bidimensionnelle de I'ycou- 
lement ytudiy 

Dans.Ie modfele proposy, les vecterus extraits fournis par le procedy d'obtention de trajectoires d'objets sont les 
45 uniques mesures utilisees pour la reconstruction du champ de vitesse dense. La ryalisation o du champ de mesure / 
= { y^, s e S} est constituye par : 



so 
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Vseao^= V^=(u^»^^) (4) 

Ainsi I'ensemble des mesures est dysigny par : 

o = {o^s G D} = { \/f ,s e D}, 

ou = (t/^, v^) est un vecteur bidimentionef dans lequel u^{resp,v^) est la composante de la vitesse suivant I'axe 
des x{resp.y^ de la vitesse foumie inltialement au site s et Dest un ensemble fini de sites. Les observations prennent 
leur valeurs dans A^bs-^obs Tensemble fini des configurations d'observations. 
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On cherche ^ construire ie champ de primitives X={Xg,s GS} sur le m§me grille S que la grille des pixels. La 
realisation e={0^s SS\={V^,s eS} du champ de primitives, Xest le champ global de vecteur Vitesse a reconstruire, et 
repr6sente la primitive au site s, i.e. le vecteur Vitesse au site s. Les primitives sont toutes de meme nature : le 
meme espace d'6tats A est associe a chacun des sites de S. L'ensemble des configurations de primitives est ^ = a®. 

Un modfele est 6labor6 afin de rendre compte de la physique du ph6nomene etudie. Ce modele est directement 
corr6l6 aux propriet6s physiques du fluide et proc^de de I'etude de la m6canique des milieux continue. 

Dans ce modele, I'energie contextuelle fait intervenir trols termes d'6nergie distincts, I'un representatif des forces 
de frottement, I'autre exprlmant les forces de presslon et le dernier repr6sentant une contralnte sur la valeur de la 
divergence de la Vitesse : 

^2(^) = ^frottemem (^) + ^pression (^) + ^divergence (^) (5) 

On conslddre une particule fluide situde en M(x,y) et une autre situde en un point trds voisin en M'(x+dx,y+dy) 
(figure 3). 

Gr§ce k Thypothfese Newtonienne des fluides incompressibles et aux lois de la m^canique des fluides, on peut 
formuler la force tangentlelle de frottement exerc6e sur la particule fluide en M par la particule fluide en M*. Ainsi, en 
assimilant I'^coulement k un 6coulement sur une surface et la particule fluide M ^ un element de cette surface, tendant 
a etre entramee par le fluide (notamment situd en M'), la force tangentlelle de frottement ffronementP^^^ ^^re exprlm6e 
comme suit : 

^frottement ~ (^) 

ou est la viscosity dynamlque du fluide, n est le vecteur unltaire normal k la surface en M et ej^ est la composante 
de V^normale a n. 

II convient d'introduire dans ce modele un terme d'6nergie rendant compte de ces forces de frottement. 

L'idee la plus naturelle est d'utiliser la fonction potentielle de laquelle derive la force f frottement A cause de la pre- 
sence de d6riv6es partielles dans I 'expression de ff^ttemenP potentiel correspondant k ffnttement ©st d*une grande 
complexite. AussI, afin de definir une 6nergie d'une fomne plus simple, on a cholsl d'exploiter la relation physique mise 
en lumlere par I'equatlon (6), k savoir quil existe une relation mutuelle d'Interactlon entre les composantes tangentielles 
de la Vitesse de deux particules voislnes. Le module de potentiel de frottement que I'on a elabore prenant en compte 
cette propri6t6 physique est d§crlt cl-apres. 

Soient s et t deux sites, n^^ le vecteur unitaire dirig6 de s vers I t^f un vecteur unltaire normal k n^, et d^f la distance 
s6parant les sites s et t On appelle ©n^^ {resp.e^) la projection du vecteur vltesse V'^du site s{resp.f) sur n^et e^^ 
la projection de Vg sur un t^^ (voir figure 4). On deflnit I'energie de frottement Uf^^^^,^ {e) par : 



dans laquelle le potentiel V'/toffemeii/,/*,^; cie la clique binaire {s^ considere est exprime par : 

^fmttement,{3,t} (^) = K \^ 9(dJ. (8) 

qui peut aussi s'^crire : 

^frottementfsj} (^) = ^j.4i p(cfj. {t^rl^s'^t))^ (9) 

Ce potentiel tradult le ph§nom6ne de glissement de deux elements de surface parall6lement I'un k I'autre. La 
variable Xj^ est un parametre ponderant la force de ce terme d'6nergie. [i est la viscosity dynamlque du fluide et p(d^^) 
ponddre I'lnfluence du terme d'6nergie en fonction de la distance des sites s et t 
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De par I'expression de I'energte de frottement, plus la viscosity est ^lev^e et plus les particules sont proches, plus 
les composantes tangentielles des vltesses de sites voisins doivent etre proches. 

Plus qu'une formule analytique precise de Up^^ssion^ on souhaite ddflnir une dnergie representative des phdnomd- 
nes de pression intervenant entre particules fluides. Ainsi, on congoit que la composante normale de la vitesse d'une 
particule ait une forte correlation avec ia composante normale de la Vitesse d'une particule situ6e sur sa normale. On 
decide d'imposer une contrainte de douceur locale forte sur les projections e^\, des vitesses de deux particules fluides 
normales (s^pardes par un vecteur colin^aire n^). On d^flnit I'energie de pression Upf^^Je) par : 

IS dans laquelle Vp^jg„ (g tj (e) est d6fini par : 
20 soit ; 

^pr9sslon,{s,t} (®) ~ 

A..H.p(rf^).(n^(V^-V,))2 (12) 

25 Des informations complementaires Ii6es aux propri6t6s physiques du ph6nomene sont disponibles et peuvent etre 

Introduites. 

Une equation fondamentate de la mecanique des fluides est la lol de conservation de la mati^re qui peut s'ecrire : 
30 |e+pV.V + Wp=0 (13) 

dans iaquelle p est la densite du flulde et Vsa Vitesse. Dans le cas ou le fluide est incompressible, c'est-^-dire que sa 
densite peut dtre considerde constante dans tout le volume de fluide et pendant tout le temps du mouvement, I'dquation 
de continuity prdcedemment enoncte prend une forme trds simple : 



35 



45 



so 
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V.U=0 (14) 



qui impose que, en tout point de I'ecoulement, la divergence de la vitesse est nulle. 
40 Si recoulement est bidimensionnel ou bien s'il est tridimensionnel avec une trds faible composante en z, il peut etre 

interessant d'exptoiter {'equation de continuity et d'en deduire une nouvelle contrainte pour I'energie contextuelle U e . 
Dans le cas des ecoulements bidimenslonnels ou semi-tridimensionnels, un temne dit de divergence Uf^fyergencei^) P®"^ 
etre ajoute k I'energie contextuelle : 



oD est un parametre ponderant influence du terme de divergence par rapport aux autres termes. 
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Revendicatlons 

1. Procddd de reconstruction d'un champ dense de caractdristiques assocides h un phdnomdne physique h partir 
de mesures ^parses. caractdris6 en ce qu'il comporte les stapes suivantes : 

une premiere 6tape (10) d'6laboratlon d'un module de reconstruction du champ elabor6 ^ partir des connais- 
sances a priori sur le phenom^ne etudi^ de fagon ^ obtenir une fonction analytique ; 
une seconde dtape (11 ) de realisation d'un maillage du champ; 

une troisi^me dtape (12) de reconstruction du champ par optimisation de la fonction analytique ; 
et en ce que le module de reconstruction prend en compte les caracteristiques physiques du ph^nomdne. 

2. Proc6d6 selon la revendlcation 1 , caract6ris6 en ce que le maillage est adapts k la repartition spatiale des mesures. 

^5 3. Precede selon la revendication 1, caracteris6 en ce que dans la seconde etape on utilise un algorithme de type 
Kohonen dans lequel afin d'assurer une bonne quality de regularisation sur le maillage, les mailles de celui-ci sont 
aussi regulidres que possibles ; et en ce que des contraintes sont imposees sur les angles entre ressorts d'une 
mallie : les quatre angles formds par les ressorts pris deux par deux devant §tre superieurs ^ un angle minimal. 

30 4. Procede selon la revendication 3, caracterise en ce qu'on impose une distance minimale inter-noeuds. 

5. Precede selon la revendication 1 , caracterise en ce que les caracteristiques sont des caracteristiques de trajectoire 
d'objet telles que la Vitesse ou I'acceieratiion. 

3S 6. Procede selon la revendication 1, caracterise en ce que les mesures eparses sont fournies par un procede de 
trajectographie d'objets. 

7. Precede selon la revendication 6, caracterise en ce que le procede de trajectographie d'objet comprend les etapes 
suivantes : 

40 

une etape d'enregistrement de signaux composes de fragments caracteristiques des positions des objets d 
differents instants, de parties qui presentent la mSme nature que les fragments mais qui sont du bruit et de 
parties attribuables au bruit sans confusion possible ; 

une etape d'extraction de parties de signaux ayant la meme nature que les fragments et determination de 
45 caracteristiques leur etant associees, ^ partir des connaissances a priori sur leur nature et sur les trajectoires 

recherchees ; 

une etape de repartition de {'ensemble des parties de signaux precedemment extraits en classes, chaque 
classe representant une trajectoire potentlelle et comprenant un nombre predetermine de fragments ; 
une etape de selection d'un sous-ensemble de classes qui satisfait des contraintes liees au type de reseau 
50 observe, comportant : 

• la generation d'un ensemble de n-uplets constitues d'un sous-ensemble de n classes, 

• la mesure de la compatibilite des classes assembiees en n-uplets par une premiere fonction analytique, 
determinee h partir des contraintes liees au type de ph6nomene observe et au dispositif d'acquisition, 

55 * . mesure de la qualite des n-uplets par une seconde fonction analytique, determinee d partir des con- 

traintes liees au phenomdne observe et au dispositif d'acquisition ; 

une etape de selection parmi les classes representant les trajectoires potentielles de celles satisfaisant les 
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contraintes. par un precede d'optimisation sous contrainte du (ou des) crit6re(s) pr6cedent(s), en utilisant les 
mesures precedentes de compatibilite et de quality, de fa^on k obtenir les trajectoires "r^elles" des objets. 

8. Precede salon Tune quelconque des revendications precedentes, caract6ris6 en ce qu'il est utilise dans le domaine 
de la velocimetrie par Imagerie de particules en nnecanique des fluides. 

9. Procedd selon la revendication 8, caract^rise en ce que la fonction analytique est telle que : 



• et designent respect ivement les mesures de Vitesse aux sites s et f du champ, 

• fg et ff designent respectivement les primitives (valeurs de Vitesse) k reconstruire aux sites s et f du champ. 

• fx est la viscosity dynamique du ftuide consider^, 

• XJ_, X\V XV et p{dgf) sont des param^tres, 

• Dest {'ensemble des mesures, 

• S est I'ensemble des sites du champ, 

• n^^ est le vecteur unitaire dirig6 de s vers f. f^est un vecteur unitaire normal k n^et ofg^est la distance s6parant 
les sites s et t, 

• VYest la divergence de la fonction f. 

10. Precede selon la revendication 9. caracterise en ce qu'une fois la fonction analytique E(/) definie, la solution f du 
probl^me est obtenue lorsque E(/) est optimisee, i.e.: 



10 




(3) 



IS 



ou : 



^ = argument optimum de f^pE(f) 



35 



ou Fest I'ensemble des solutions possibles. 
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DONNEES EPARSES FOURNIES PAR 
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RESOLUTION DU PROBLEME D'OPTIMISATION 
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